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The peculiar correlation effects on the triangular lattice are studied by means of the rotation- 
ally invariant slave boson method in a cellular cluster approach. Hence nonlocal correlations are 
included in a short-range regime. Their impact for the single-band Hubbard model is studied at 
half filling, i.e., on the Mott transition, and with doping. Using the realistic band structure of 
Nai:Co02, we may also shed light on the cobaltate physics for 2;>l/3, with the in-plane transition 
from antiferromagnetic tendencies towards the onset of ferromagnetism for a finite Hubbard U . 
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The transition- metal oxide Naa;Co02, consisting of 
C0O2 layers and Na ions inbetween, is one of the most 
fascinating examples of a doped correlated electron sys- 
tem on the triangular lattice. Albeit evidently the bare 
C0O2 (a;=0) compound is metallic surprisingly, the 
effect of strong correlation appears to be severe close to 
the band-insulating (x=l) limit. Several experimental 
studies have revealed a rather rich phase diagram with 
highlighted physical properties such as superconductiv- 
ity close to a;=0.3 (when intercalated with H2O) 
Pauli(Curie-Weiss)-like mctallicity for x<0.5(>0.5) 
in-plane antiferromagnetic (AFM) order at a;=0.5 
large thcrmopowcr around 0.71<a;<0.84 charge dis- 
proportionation for a;>0.5 @ and in-plane ferromagnetic 
(FM) order for 0.75<a;<0.9 0, Q- Due to the sizable 
(^2g,eg) crystal-field (CF) splitting, the Co ion is ex- 
pected to be in a low-spin state, with x controlling the 
residual occupation of the t2g manifold. Hence Co^+ 
{S=l/2) for x=0 and Go^+ (S'=0) for x=l. Yet the 
respective t2g fillings and apparent Fermi surface (FS) 
for smaller x arc still a matter of debate 0]. Though 
calculations based on the local density approximation 
(LDA) yield an f2g-iiiternal aig-e'g CF splitting of -0.1 
eV [10|, lul, only the aig-like bands are expected to be 
partially depleted for larger x. Those should form a sin- 
gle hexagonal FS sheet that is hole- like, i.e., the hop- 
ping t within a nearest-neighbor (NN) tight-binding (TB) 
model is negative. 

Concerning magnetism, LDA suggests FM order for al- 
ready small X, although the AFM state is very close in 
energy Dynamical mean-field theory (DMFT) 

studies for the NN-TB Hubbard model on a triangu- 
lar lattice [13, show an instability towards FM or- 
der only for t>0. Gao et al. described the appear- 
ance of a renormalized Stoner instability in an infinite- 
U Gutzwiller treatment of a third NN-TB model for 
Na^Co02 (with the NN t<0) at X'^0.67. Recently, a 
finite-?/ LDA-I- Gutzwiller approach obtained intralaycr 
FM order at larger x , but contrary to experiment [l'] 
finds full FM order at low doping. There arc cluster 
approaches to the Hubbard model on the triangular lat- 



tice (e.g. [iilii) but without incorporating the detailed 
electronic structure of Na2,Co02. 

In the present work, realistic sodium cobaltate is in- 
vestigated at larger Na doping, finite U and by explicitly 
including NN correlation effects via the cellular cluster 
scheme. The puzzling change from AFM tendencies at 
small X towards the onset of intralayer FM order may 
be described within a Hubbard-like model using LDA 
dispersions and moderate U. The recently generalized 
rotationally invariant slave boson mean-field theory [2^ 
(RISB) is applied to the problem [sO], thereby tiHng the 
lattice into NN triangles (see Fig. [1]). To benchmark this 
approach in the present setup, we first studied the NN- 
TB single-band Hubbard model, written in the cellular 
cluster scheme as H=H(K) + Ha, i.e.. 



H 



aijcr 



Tin 



(1) 

where a marks the cluster, ij are site indices on the clus- 
ter, cr denotes the spin (f, i) and K is the cluster wave 
vector. Note that for the cluster dispersion J^K^^acr—^ 
holds and the intraclustcr hopping is taken care of by 
the second term in ([T]). Within RISB the electron op- 
erator CiCT is represented as Ci^—R[(j)]f-^ fja', where R 
is a non-diagonal transformation operator that relates 
the physical operator to the quasiparticlc (QP) opera- 
tor fi^. The transformation R is written in terms of 
slave bosons {(fiAn} with two indices, namely A for the 
physical-electron state and n for the QP Fock state. It 
follows that the kinetic part of ([1]) is expressed via the QP 
operators with renormalized dispersions and the operator 
character of the local part is carried solely by the slave 
bosons [2^. Two constraints, namely J2 A71 '^lin'f' An~'^ 

and J2Ann' An''t^A,M\flfja'W)^fl fja'^ tO SClect the 

physical states are imposed through a set of Lagrange 
multipliers {A}. A saddle-point solution is obtained by 
condensing the bosons and extremalizing the correspond- 
ing free energy F{{(f)}, {A}) [3l[. The nonlocal QP weight 
matrix at saddle point reads Z=RRt. 

In the following, we study the paramagnetic (PM) 
phase where x denotes the electron doping normalized 
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FIG. 1: (color online) 3-site cluster approach on the triangular 
lattice for isotropic NN hopping t. Left: lattice tiling, right: 
QP DOS for different interaction strength and doping. 



to a single orbital. The QP density of states (DOS) 
of the noninteracting model (bandwidth W=9t), along 
with comparative interacting cases, is shown in Fig. [1] 
As displayed in Fig. ^ the RISB treatment for x~0 re- 
veals a first-order Mott transition at Uc^l2.2t=1.36W , 
in agreement with values from exact diagonalization plj 
(-12.lt), cellular DMFT ^ (-10. 5t) and variational 
cluster approximation [Hi (~12i). Upon doping, the 
metallic state may be extended to larger U, although a 
breakdown of the conventional PM phase is found for 
small X at U>Uc- The magnitude of the offdiagonal 
QP weights Zij as well as of the NN spin correlation 
(SiSj) are enhanced close to the Mott transition. Note 
that (SiSj) is negative, i.e., of AFM character. Close 
to half filling, the Zij depend strongly on x and the on- 
site spin correlation is substantially enhanced for U^O. 
The AFM strength of the intersite (SiSj) is decreasing 
with a;, but remains slightly above the U=0 case in the 
(meta)stability range of the PM phase. Note that the 
global FM instability was located around [3]. 

For a realistic description of Naa;Co02 we first per- 
formed LDA calculations at different stoichiometrics. A 
supercell involving the NN Co triangle of a given C0O2 
layer served as the base structure. Decorating the lat- 
ter with additional Na ions above the Co sites yields 
the dopings a;=0, 1/3, 2/3,1. No sodium ions above/below 
oxygen positions are considered, which should be ade- 
quate for the overall behavior with doping. A minimal 
model for the band structure is derived by represent- 
ing the low-energy states by three maximally-localized 
Wannier functions [2§| (WFs) of aig type (see Fig. 
Thus the full dispersions (M^~'1.2 eV) are approximated 
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like 



through a downfolding procedure [23| 
bands. Allowing for the essential behavior with x, a dop- 
ing dependent hamiltonian is obtained through a linear 
interpolation HldaIK, x)^xHl^j^{K)+{l-x)Hl^j^(K), 
where Xa,Xb are the neighboring LDA-treated dopings. 
Note that the K points are generally on a 3-dimensional 
(3D) mesh, describing the full band dispersions. In the 
full hamiltonian ([T]) the kinetic term is now given by 
-ffLDA(K,x) without its onsite (cluster) term t^. The 
latter replaces the model t in the "onsite" quadratic term 



of ([T]), incorporating both, the intersite hoppings on the 
Co cluster and the CF-terms on the individual ions. Be 
aware that the level of frustration is diminished via the 
symmetry-breaking Na occupations for intermediate x. 

Contrary to the former model study, the realistic NN 
hoppings are now negative 1J|. Their absolute value is 



more or less linearly decreasing with x, in average from 
0.21 eV (a;=0) to 0.11 eV (a;=l). Any further explicit 
exchange couplings J are neglected. Integrating out the 
O states should render such couplings necessary (e.g., 
(A)FM exchange between the Co ions). Yet a rigorous 
computation for metals is hard [isj and experiments re- 
vealed Q that the values for a pure spin hamiltonian are 
small. Calculations including Heisenberg-likc terms in 
H resulted in no qualitative change for reasonable val- 
ues of J (see footnote on the spin susceptibility). No 
long-range Coulomb potential is included, hence explicit 
charge-disproportionated states do not appear (though 
site occupations surely differ). We concentrate the dis- 
cussion on the phases for a;>l/3, since neglecting the 
remaining e'g orbitals at small x may not be improper. 
In fact, when performing the RISB calculation at a;=0, 
the Mott transition sets in at ?7c~ 3.1 eV. But this value 
is surely too small because of our orbital restrictions. 

Figure m shows RISB results for a;=l/3,2/3 and with 
varying Na content for reasonable ranges/values of the 
interaction strength U. The resulting moderate size of 
the diagonal QP weight is in good agreement with exper- 
iment. Note that for a;=l/3 the offdiagonal Zy between 
sites on the triangle is positive with a maximum at U'^S 
eV. For a;=2/3 its value is negative and still increasing in 
magnitude in the studied U range. The sign change im- 
plies that the hierachy in the correlation strength within 
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FIG. 2: (color online) Observables for the NN-TB Hubbard 
model {t>0) on the triangular lattice. Left: U dependence 
for (a) QP weights, (b) onsite spin correlation and (c) NN 
spin correlation. Right: doping dependence of the same ob- 
servables. The Zij are multiplied by 5 for better visibility. 
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FIG. 3: (color online) From left to right: aig-Iike WFs in 
Naa;Co02 unit cell for 1=2/3 on Co(l), Co(2) and Co(3). 
Large blue (dark) balls denote Na ions, small red (gray) ball 
mark O ions. 
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FIG. 4: (color online) Observables for Naa;Co02. Left: (a) 
site-resolved QP weights, (b) onsite- and (c) NN spin corre- 
lations for fixed doping. a:=l/3: (black/dark) lines, a;=2/3: 
(red/gray) lines. The offdiagonal Zij go to zero at ?7=0. Co 
sites are numbered with Na filling: '1' has Na on top at a::=l/3, 
'2' also at a;=2/3. Right: doping dependence of the averaged 
observables. All Zij are multiplied by 50 for better visibility. 



the bonding-antibonding cluster orbitals is reversed. The 
onsite spin correlations are saturating for large IJ , while 
the intersite (SiSj) are always negative but decrease in 
magnitude in the latter regime. For a;=l/3 the NN (spin) 
correlations are enhanced in magnitude between Co sites 
without Na on top. Due to CF effects, the orbital fill- 
ing on a Co site tends generally to increase when placing 
Na on top. For the doping study, ?7=3 eV and U =5 eV 
were used. There the intersite Zij changes sign at x<0.5 
and has a minimum close to a::=2/3. As expected, the 
local spin correlations are strongly increased for finite U 
close to half filling The AFM character of the NN 
spin correlations is strengthened for intermediate x com- 
pared to the [/=0 case. However interestingly, there is a 
crossover at x~2/3 where the AFM magnitude for C/t^O 
becomes lower than in the latter case. Thus the AFM 
tendencies are surpressed by electronic correlations for 
large x, where the superexchange is diminished. This 
seems reminiscent of the Nagaoka mechanism 2J1 for the 



infinite-[/ Hubbard model on the square lattice, leading 
to FM order. But the Nagaoka state is not a g round state 
on the triangular lattice for a single hole [25j . 

To elucidate the problem of magnetic instabilities, a 
residual magnetic field Hf^J/ie^ is applied in the PM 
phase. Therefrom the uniform magnetic susceptibility 
is determined as x=9A//9-fff«(5Af /(5/i, where SM is the 
residual magnetization induced by 8h. Figure [5^ shows 
the normalized response x/xoi where is the suscep- 
tibility for t/=0 [34]. It is seen that there is strong re- 
sponse for a;~2/3 and C/=5 eV, similar to results in the 
infinite- [/ limit Additionally, there is already a pre- 
cursive regime for a;>0.62 where FM ordering tendencies 
show up. This coincides with recent NMR measurements 
by Lang et al, who find a crossover from AFM to FM 
correlations at a;=0.63-0.65 [2^. From the site- and com- 
ponent resolved intersite spin correlations in the applied 
field, the dominant response of {SizSjz) in the doping 
regime 0.62<a;<0.7 is again obvious. Note that also J7=3 
eV exhibits minor FM ordering tendencies there. Inter- 
estingly, for a;~0.35 one may observe intralayer AFM re- 
sponse for both values of U . An AFM ordering signal 
is identified through locally favorable spin antialignment 
between Co(l)-Co(2) and Co(l)-Co(3), whereas Co(2)- 
Co(3) favor spin alignment. Remember that for a;=l/3 
the Na ion is above Co(l), thus the Co differentiation 
via neighboring Na ions obviously triggers the AFM ten- 
dencies. The outer-field induced AFM response is in line 
with the observation that such a field can lift the ef- 
fects of kinetic-energy frustration on a triangular clus- 
ter [S^l ■ The fact that no clear long-range order is visible 
for larger x (i.e., negative x/xo) may be explained by 
the present 3D model that cannot stabilize the A-type 
AFM phase (no proper interlayer resolution). Yet it is 
possible to stabilize FM order within a single layer by 
using only the 2D-projected LDA dispersion. Figure [5f 
shows the tie-line construction between the PM and FM 
phase in this case. A stable homogeneous FM order is 
revealed for x>0.74, with magnetic moments close to ex- 
perimental values. The onset of the heterogeneous phase 
at a;'-^0.6 matches the begin of FM tendencies in the 
full 3D model. In this context, hints for a spin-liquid 
ground state at 0.71<a;<0.75 [28| are rather interesting. 
Either 2D phase becomes hard to stabilize numerically 
for x>0.87. Though at a;=2/3 the Fermi level is located 
in the upper maximum peak of the QP DOS (see Fig.|4l5), 
our study reveals that a Stoner instability can not be the 
sole origin for magnetic ordering. The nonmonotonic be- 
havior of Zij, the subtle change of the frustration level 
with X, the precursive FM regime, the PM-FM phase 
competition as well as the breakdown of the 2D phases 
at very large doping are all pointing towards nonlocal 
correlations as an additional important ingredient. As 
a side effect, plotting d'^F/dx'^ (Fig. S^-) shows that the 
present 3D-PM phase becomes unstable with respect to 
doping for x>0.7, while U stabihzes the phase for a;<0.5. 
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FIG. 5: (color online) Ordering tendencies and phase stability 
in Naa:Co02. (a) normalized magnetic susceptibility, (b-d) 
component-resolved NN spin correlation in the small field Sh. 
Solid black line: (Si^Sj^)sh for U=0. (e) QP DOS for x=2/3. 
(f ) RISB free energies for PM and FM phase of the single layer 
for U—5 eV. Straight cyan (lightgray) line: common tangent. 
Inset: net magnetic moment per Co ion in the FM phase, (g) 
Second derivative of the full PM free energy. 

Phase separation is experimentally indeed observed at 
large x (e.g. [2^ and references therein). 

In conclusion, it was shown that the RISB formalism 
may describe the essential physics of the NN-TB Hub- 
bard model on the triangular lattice, including the first- 
order Mott transition and the Fcrmi-liquid upon doping. 
The magnetic behavior with doping in Na2;Co02 can be 
understood with a cluster Hubbard model using LDA 
dispersions. In very good quantitative agreement with 
experimental data, the change from AFM to FM ten- 
dencies with a final first-order transition in x into an in- 
plane FM phase was revealed. For the correlated physics 
with intralayer FM tendencies starting at a;>0.62 and 2D- 
FM order for 0.74<a:<0.87, an interaction strength C/>3 
eV is sufficient. Since then U/W>2.5 and the strongly- 
correlated regime is already reached, without implying 
that x=0 is Mott insulating (due to the probably in- 
creased orbital fluctuations there). Further studies, also 
at small doping, including long-range and inter layer cor- 
relations are needed to clarify more details. 

The author is indebted to M. Potthoff, I. I. Mazin and 
P. S. Cornaglia for helpful discussions. This work was 
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